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Working with Lieb's transfer matrix for the dimer model, we point out that the full set 
of dimer configurations may be partitioned into disjoint subsets (sectors) closed under the 
action of the transfer matrix. These sectors are labelled by an integer or half-integer quantum 
number we call the variation index. In the continuum scaling limit, each sector gives rise to a 
representation of the Virasoro algebra. We determine the corresponding conformal partition 
functions and their finitizations, and observe an intriguing link to the Ramond and Neveu- 
Schwarz sectors of the critical dense polymer model as described by a conformal field theory 
with central charge c — ~2. 
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I. INTRODUCTION 



In his 1967 paper 



II, Lieb proposed a transfer matrix for the dimer model 



m 



It is our impression 



that this work has not been fully exploited and that it may shed light on the status of the dimer model 
as a conformal field theory (CFT). A certain confusion has prevailed for some time regarding the central 



charge of the relevant CFT. Indeed, some aspects of the dimer model can be described 
with c = -2 j?- 



12l |. hintin g at its logarithmic character [I3l4l6l|. while some other aspects, most notably 



S] by a CFT 



those involving monomers jl7-19], seem to require 



21[ a description in terms of a CFT with c = 1. 



Moreover, in certain geometries, the dimer model is equivalent 



22 



23l | to the Abelian sandpile model 



26l |. which is widely believed 27h31I] to be described by a logarithmic CFT with c = —2. 

Here we consider dimers defined on rectangular lattices with or without periodic boundary conditions 
on the horizontal or vertical edges. Working with Lieb's transfer matrix, we find that the full set of dimer 
configurations may be partitioned into disjoint subsets we call sectors. These sectors are closed under the 
action of the transfer matrix and are naturally labelled by an integer or half-integer quantum number we 
call the variation index. In Lieb's spin picture, this quantum number is obtained as the eigenvalue of a 
particular matrix defined in terms of Pauli matrices. 

In the case with open boundary conditions, each sector is found to be associated with a representation 



* Electronic address: 'j. rasmussen@uq. edu. au| 
^Electronic address: ,pliilippe.ruelle@uclouvain.be| 



2 



of the Virasoro algebra as the corresponding sector partition function yields a Virasoro character in the 
continuum scaling limit. These characters can be organised in a manner consistent with a CFT with 
central charge c = —2. Furthermore, the corresponding finitizations are accounted for by a combinatorial 
enumeration reminiscent of the physical combinatorics underlying the critical dense polymer model on 



the strip 
Following 



11 



and cylinder 



12l |. This gives rise to a separation into Ramond and Neveu-Schwarz sectors. 



32l | , we are thus using the terminology of supersymmetry even though we do not claim any 



superconformal symmetry in the dimer model. 

In the case of even system sizes with periodic boundary conditions, each sector can be described 
as composed of two 'chiral halves', much akin to the critical dense polymer model on the cylinder [12 1 
and many other similar lattice models. The two 'chiral halves' are not independent, but matched up 
by imposing a simple gluing condition. The combinatorial description of the ensuing finitized partition 
functions again resemble the physical combinatorics of the critical dense polymer model. 

Despite the obvious parallels between our treatment of the conformal properties of the dimer model 
and the similar studies of the critical dense polymer model 
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VM, there are crucial differences. In 



particular, the variation index is a good quantum number, while the related so-called defect number 
in [ill, 3] is not. In addition, the variation index takes on negative as well as positive values, again 
unlike the defect number which is non-negative. These differences ensure that the physical partition 
function for the dimer model on a torus with both lattice sides even is modular invariant, as opposed to 
the toroidal partition functions discussed in [12]. This modular invariant partition function for dimers 
was originally obtained by Ferdinand [s^ , and has since appeared as the partition function for symplectic 
fermions as described by the so-called triplet model [lo| . 

The layout of this paper is as follows. In Section [Til we review the basic formulation of the transfer 
matrix, as proposed by Lieb [l[. Section IlIII examines the partitioning of the configuration space into 
an (in the continuum scaling limit) infinite number of sectors left invariant by the transfer matrix T. 
These sectors correspond to orbits under T and are thus the finest possible partitioning into invariant 
spaces. The sectors are interpreted in terms of eigenspaces of a certain operator V, whose main property 
is to commute with the squared transfer matrix T^. Section ITVl generalizes Lieb's diagonalisation of the 
(squared) transfer matrix to the open boundary conditions, and works out the relevant partition functions. 
Sections [V] and [Vlj extend this to the individual sectors, whose conformal contents are discussed. This is 
made more explicit in Section [VIII and IVIIIl respectively for the open and periodic boundary conditions, 
where the combinatorics of the conformal states is related to that in the critical dense polymer model 
Concluding remarks are presented in Section llXi Jacobi's theta functions are reviewed in 
Appendix [Aj while Appendix [Bj contains some technical details. 
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Note added. 

After completion of the present work, Jesper L. Jacobsen has informed us that the existence of sectors 
in the dimer model had been observed before, albeit not in the framework of Lieb's transfer matrix, nor 
applied to the computation of conformal partition functions. Two descriptions have been proposed. In 
the first one, the separation into sectors is based on the geometric observation that the superposition 
of a general dimer configuration and a certain fixed configuration results in the formation of a number 
of 'strings'. In this view, a sector collects all dimer configurations giving rise to the same number of 



strings. The details of this construction can be found in 
has been discussed in the context of trimer tilings 



35|. The other, more algebraic approach 



3g], but works more generally for p-mers, p > 2. In 



this second description, a certain function on p letters is associated with a 'profile', namely a line that 
goes through the lattice from left to right without crossing any of the p-mers. This function is constant 
on the set of all profiles which can be obtained from each other by adding p-mers (i.e. the portion of the 
lattice in-between two profiles can be covered by p-mers) and defines the sectors. For p = 2, this function 
essentially coincides with our variation index. Despite these earlier observations, we believe that our 
applications of Lieb's transfer matrix provide new insight into the sector structure as well as the dimer 
model more generally. 



II. DIMERS AND LIEB'S TRANSFER MATRIX 
A. Dimer model 

Let us start by briefiy recalling what the dimer model is. We consider a rectangular grid in I? with 
M rows and N columns (possibly with periodic boundary conditions in the vertical and/or horizontal 
directions), which we cover with dimers. A dimer is a rectangular tile (1 x 2 or 2 x 1), like a domino, 
covering exactly two adjacent sites^. A dimer configuration is then an arrangement of dimers such that 
every site of the grid is covered by exactly one tile, thereby disallowing dimers to overlap. A primary 
objective in this model is to study the statistics of the dimer configurations. 

The configurations are weighted according to their tile content. Without loss of generality, we may 
decide that a vertical dimer has weight 1, while a horizontal dimer has weight a, so that a generic 
configuration receives a weight equal to a^, with h the number of horizontal dimers. The partition 
function sums the weights of all configurations, 

Zm,n= Y1 (2.1) 

configs 



There is a well-known generalisation in which monomers, which cover only one site each, are allowed, but we do not 
consider this case here. 
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B. Transfer matrix 

In Lieb introduced a transfer matrix for this model, which builds all possible dimer configurations 
row by row with the direction of transfer upward. First, one replaces a dimer configuration by a set of 
arrows, up or down, attached to the sites, and placed above them, as illustrated in Figure 1. An up 
arrow at site i means that there is a vertical dimer covering i and its northern neighbour i + e^; a down 
arrow indicates the absence of such a vertical dimer. In the latter case, it either means, if the down 
arrow is right above an up arrow, that site i pairs up with its neighbour below to form a vertical dimer; 
otherwise, it means that the down arrow at i associates itself with another down arrow at the left or right 
neighbouring site, to form a horizontal dimer. This correspondence is indicated in Figure 1 for a tile 
configuration covering five successive rows. A row of N sites can have 2^ different arrow configurations. 
It is natural to think of the up and down arrows as the canonical base elements (1,0) and (0,1) of a 
vector space C^. A row configuration of arrows can then be viewed as an element of the multiple tensor 
product (C^)®^. 

The next task is to define the operators which produce dimers; vertical or horizontal. Their action 
will be given in terms of tensor products of Pauli matrices, 

o-i = I(g)...(g)][®fT(g)I(g)...(g)I, (2.2) 

where a, in position i, is one of the following three matrices (given in the basis |t), |i)) 

/O 1\ 



1 





■ • 


/ 










:) 




0/ 


Vo 





(2.3) 



The matrices on different sites commute, fjjfjj = ajai for all i ^ j, while those referring to the same site 
satisfy (<t~)^ = (it+)^ = (fif )^ - I = and af af of = af and hence af af af = af. 
The transfer matrix acts on a horizontal sequence of arrows in two steps. 




FIG. 1: An arrow configuration is associated with a dimer configuration: a vertical dimer at a site i becomes an 
up arrow (in brown) placed above i, while the absence of a vertical dimer is indicated by a down arrow (in grey). 
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1. It first reverses all incoming arrows, by acting with 

= n (2-4) 

i 

This ensures that an up arrow becomes a down arrow, in accordance with the fact that an up 
arrow cannot propagate since two vertical dimers would then overlap. Similarly, a down arrow 
automatically propagates upward to an up arrow, corresponding to the lower half of a vertical 
dimer. 

2. After step 1, an outgoing up arrow at site i, corresponding to a vertical dimer, may be kept or 
changed into half a horizontal dimer. A horizontal dimer can be produced by reversing a pair of 
neighbouring up arrows, i.e. by acting with cj^a^j^^^ multiplied by the corresponding weight a. The 
combination I + aa^a^j^-^ does not do anything or implements the change, so the operator 

^3 = n P + " ^^i'^i+i] = exp (a c^ioi+i) (2.5) 

i i 

offers this alternative for every pair of adjacent sites in a row. 
The product V3V1 incorporates all possible ways to go from a row of specified arrows to the next one 



rerever possible, and correctly handles 
2 



up while respecting the constraints: it creates horizontal dimers w. 
the incoming vertical dimers. It is the transfer matrix defined in 

TV 

T = V^V, = exp {J2 « ^r^r+i) n • (2.6) 

i>l j=l 

If no horizontal periodicity is assumed, the summation over i in the exponential runs from 1 to — 1, 
whereas in the periodic case, it runs to N, with the identification ciAr+i = cxi. 

The action of T on an incoming row of arrows |in) produces a linear combination of outgoing rows 
of arrows, the set of which represents all possible row configurations atop |in). The coefficients in this 
linear combination involve powers of a which count the number of horizontal dimers. Likewise, the power 
T"^ iteratively constructs all possible arrays of m rows built on top of |in). The top row in this process, 
|out) = r^lin), comprises all possible outgoing row configurations, m layers higher than the incoming 
row I in). The configuration in Figure 1 represents a possible realisation of the action of on the bottom 
row, taken as incoming state |in) = | tiittiiitt) • I^i this particular case (A = 10 with no horizontal 
periodicity), there are 24 203 other possible realisations compatible with |in). 

We note that, although the mapping of a dimer configuration to an arrow configuration is well-defined 
and unambiguous, the inverse mapping is locally not so. The source of ambiguity resides in the pairs 



^ Lieb's original matrix also allows for monomers in which case a third operator is needed in addition to Vi and V3. In his 
paper [J], this operator is denoted by V2. 
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of adjacent down arrows, which, locally, may correspond to a horizontal dimer or to two vertical dimers 
coming from below. The ambiguity can, however, be lifted by looking at the neighbourhood of these 
pairs: a pair of down arrows atop a pair of up arrows unambiguously corresponds to two vertical dimers. 
Thus, by looking at larger and larger neighbourhoods, and eventually at the whole grid, the mapping 
from arrows to dimers becomes globally unambiguous for a prescribed set of boundary conditions. 

We further emphasise the important role played by the boundary conditions as they may forbid certain 
rows of arrows. For example, the bottom row of arrows in Figure 1 would not be allowed if that row had 
been a straight boundary: the three adjacent down arrows in the right half necessarily give rise to an 
arrangement of dimers where a dimer protrudes from the grid. 

C. Partition functions 

The partition functions for various boundary conditions can be computed in terms of T in the usual 
way. If vertical periodicity is imposed, we have 

Zm,n = TtT^' (2.7) 

for periodic or non-periodic (open) horizontal boundary conditions. In case no vertical periodicity is 
imposed, the grid is bordered by two horizontal edges. On these edges, the usual condition is to restrict 
the dimers not to protrude from the boundaries. In this situation, the partition function is given by 

ZM,N = Y.{ont\T^'-'\m), (2.8) 

|in> 

where the sum is over all arrow configurations |in) allowed on the lower boundary, namely those where 
down arrows come in adjacent pairs. The out-state is given by 

|out) = III • • • I) (2.9) 

since the top row cannot contain vertical dimers sticking out of the grid. The partition function ()2.8p is 
thus given by the sum of certain elements of the row corresponding to the arrow configuration |out) in 
the matrix T^'^~^ . An obvious advantage of using the transfer matrix is that it allows, through formulae 
such as (|2.8p . to consider other, new boundary conditions on the bottom and top edges by choosing other 
sets of in- and out-states, thereby allowing the dimers to protrude in a specified way. 

By construction of the pure dimer model, the partition function Zm,n vanishes if both M and N are 
odd. For convenience, we shall therefore assume M even in the following. As we are mainly interested in 
situations with vertical periodicity as in (12. 7p . our analysis is primarily based on instead of T itself. 
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III. SECTORS AND VARIATION INDEX 



The transfer matrix is typically not irreducible on the space of arrow configurations; it has many 
invariant subspaces. If such a subspace is spanned by a set of arrow configurations, we refer to it as a 
T-sector. Note that a generic eigenspacc of T is not a T-scctor, since it is spanned by a set of linear 
combinations of arrow configurations where the number of participating arrow configurations (generally) 
exceeds the dimension of the eigenspace. Similarly, a T^-invariant subspace spanned by a set of arrow 
configurations is called a T^-sector. 

The simplest example of a T-sector is the subspace spanned by the two alternating sequences, 

imt---) and Itm---), (3.1) 

where must be even in the case of periodic boundary conditions. (Periodic boundary conditions with N 
odd are discussed below.) As the transfer matrix maps one of these two states onto the other, they clearly 
form a two-dimensional invariant subspace. Individually, the two configurations form one-dimensional 
T^-sectors. The other sectors are larger and have dimensions which grow with N. 

In the following, we will determine the finest possible decomposition of the set of arrow configurations 
into T- or T^-sectors, and work out the corresponding dimensions. This classification depends on the 
parity of N and whether periodicity is imposed. 



A. Variation index 



Let us introduce the operator 



f N 



N 



N 



V = E (-1)' k+ar - 1] = 1 ^(-1)V,^ = 



i=l 



i=l 



^ (- 1)* ajar if N is even, 

i=l 

N 

i + ^(-l)V+ar if AT is odd. 



(3.2) 



i=l 



Each of the constituent operators 

(-1)' \atar -\]= ^^al 



where a^ 



1 
-1 



(3.3) 



takes the value +| or —\ depending on the parity of i and on whether the arrow at i is up or down. V 
is thus diagonal in the arrow basis and has eigenvalues v of the form 

A^ 



V G {— n, — n -1-1, . . . , n}. 



n 



(3.4) 



We note that a pair of equal neighbouring arrows in positions i and l<z<A?^, brings a zero 

contribution to while two opposite neighbouring arrows in similar positions contribute or —1 to v. 
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We may therefore interpret v as an accumulative measure of the local changes in an arrow configuration; 
for this reason, we call it the variation index. This variation index v is an integer for N even, but a 
half-integer for odd. 

In the case of periodic boundary conditions, the variation index depends on where the labelling i of 
the arrows begins. For N even, an odd shift in the labelling will change the sign, whereas an even shift 
has no effect on the variation index. For N odd, on the other hand, the starting point significantly affects 
the value of the variation index. 

For each value of v in ()3.4p . let Ey be the vector (sub-)space spanned by the arrow configurations 
having variation index equal to v. For instance. En is generated by the single state | itit " " ") discussed 
above (13. ip . while E-n is generated by its flip. It follows that dim En = d\mE-n = 1. More generally, 
we have 

dim£^„ = dim£^_„ , (3.5) 

since a global flip of arrows bijectively maps one of these spaces onto the other. Because each site 
contributes it^ to u , the generating function for the dimensions of E^ is given by 

(dimS,)r = Ut + ^Y , (3.6) 



where the summation index v increments in steps of 1, from which it follows that 



dim^^ = ( ^ ). (3.7) 



n — V 



B. Sectors 



Our classification of the finest possible sector decompositions of the space of arrow configurations is 
conveniently described in terms of the variation index v and the spaces E^. First, we show that certain 
simple direct sums of these spaces constitute T- and T^-invariant spaces. We subsequently argue that 
these sectors are indeed the smallest possible ones. Our analysis depends on the parity of and whether 
periodic boundary conditions are imposed. 

We first examine how the transfer matrix T = exp (^ • a(j^a^_^^) Hj ^^^^ °^ Because Vi = 
W - crj induces a global flip, it anti-commutes with V. Moreover, each constituent operator (t~^^ 
commutes with V unless i = N, N is odd and periodic conditions are imposed horizontally. Except 
in this particular case, V and T therefore anti- commute. It follows that T changes the sign^ of the 
eigenvalues of V, thus mapping E^ onto E—y, 

T: Ey ^ E^y. (3.8) 



^ Likewise, V changes the sign of the eigenvalues of T, which therefore come in pairs of opposite sign except in the sector 
Eq on which V vanishes. 
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Consequently, the following spaces are left invariant: 

T-sectors : Eq, E^®E_^, v > Q. (3.9) 

It also follows that commutes with V and thus leaves each individual space E^ invariant, 

T^-sectors : E^. (3.10) 

The T- orbit (or T'^- orbit) of a given arrow configuration is introduced as the linear span of the set 
of arrow configurations appearing in the expansions of repeated applications of T (or T^) on the given 
arrow configuration. From the explicit form of T, an element in the T-orbit of C is a linear combination 
of elements of the form 



c 



where each factor in square brackets represents the action of one T. Within a given such factor, the 
position labels of the o"~ matrices are all different. 

We note that being in a given orbit is an equivalence relation. To see this, we observe that if the 
previous relation ()3.1ip involves q factors of T, namely C = [Tq] . . . [T2] [Ti] C 7^ 0, then we also have 
C = [Ti] [T2] . . . [Tg] C because each factor squares to 1 (unless it annihilates the configuration it is acting 
on). This shows that the relation is symmetric, from which transitivity readily follows. 

As the orbits provide the maximal decomposition into invariant subspaces, it is enough to show that 
the sectors Eq and E^ © E^y in (|3.9p coincide with orbits under T. As indicated, the above proof that 
V anti-commutes with T breaks down in the periodic case with N odd. Indeed, we show below that in 
this case, there is a single orbit so that T is transitive on the full space. We start with the open case as 
it simplifies the discussion of the periodic case. 

To show that every T-sector in (j3.9p is a T-orbit, it suffices to demonstrate that a particular canonical 
arrow configuration in the sector is an element of the orbit of any arrow configuration in the sector. Such 
a canonical arrow configuration is conveniently chosen as 



Co = |ii---i) G^o, (or |tt---t)G^oj (3.12) 

and, for > 0, as 

Cv = I int - --m i--i) e Ey or = I titi - --m t--t) g e^v (3.13) 

2[v] N-2[v] 2[v\ N-2[v\ 

We now outline an algorithm for reaching one of the canonical configurations by repeated applications 
of T on a general arrow configuration. 

A general arrow configuration can be partitioned into sections consisting of aligned arrows in such 
a way that two neighbouring sections contain opposite arrows. If the configuration is composed of one 
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section only, it is already in canonical form ( (j3.12p for N even; ()3.13p for A*" odd), so let us assume it 
consists of several sections. From the identities 

= (yicyi+iViVi , (yt^yt+i = Viarar^^Vi , (3.14) 

we see that flipping a pair of adjacent aligned arrows produces another configuration in the same T-orbit 
(in fact r^-orbit), and similarly for an even length sequence of aligned arrows. If the leftmost section 
in our general configuration contains an even number of arrows, we flip these arrows. The result is a 
configuration whose (new) leftmost section is enlarged. If this section is even, we repeat the previous 
flipping step. This is repeated until we have produced a canonical configuration with all arrows aligned 
or a configuration whose leftmost section contains an odd number of arrows followed by at least one other 
section. In the latter case, we may now 'move' the interface between the two leftmost sections all the way 
to the left by application of YYj=2'^J YYj=2^~j ^ where m is the odd number of arrows in the leftmost 
section. The result is a configuration in which the two leftmost arrows are opposite. We then repeat 
the procedure just outlined on the remaining N — 2 arrows. After a certain number of repetitions, the 
resulting configuration will be of the form 

l(it)(it) • • • (it) tt • • • t) or |(ti)(ti) • • • (ti) ii • • • i). (3.15) 

If is even, the number of non-embraced arrows to the right is even, and flipping these produces a 
configuration in canonical form (j3.13|) . If N is odd, there is an even number of aligned arrows in the 
rightmost section, and flipping them produces a canonical configuration of the form (|3.13p . This completes 
the proof in the open case. 

Let us now turn to the periodic case. Because the periodic transfer matrix TP*^"" involves an extra factor 
compared to the open matrix r°P^°, namely T^^'^ = [I + a cjjyCJ]^] T°p'^", it follows that the TP'^^'-orbits 
must be unions of the r°P'^'^-orbits, that is, unions of Eq and © E^^ for appropriate values of u > 0. 

For N even, the anti-commutation of the transfer matrix with V implies that T cannot connect two 
sectors with different values of 1^;!. Therefore, the T-orbits in the periodic case coincide with those of the 
open case, and are given by the sectors (13. 9p . 

For N odd, the previous argument does not apply since the anti-commutation of T and V is lost. It 
is, however, not difficult to see that the full configuration space forms a single orbit under T. Indeed, let 
us consider the (unique) configuration |(ti) ■ ■ ■ (t4')t) iii ^-n- Since its first and last arrows are up, by 
applying cr^cr^ = a^a^ViVi, we obtain a non-zero configuration of E^n+2, itself related, by some power 
of T^, to another configuration of E^n+2 with its first and last arrows up (if —n + 2 < n — 1). Iterating 
these transformations, we see that the subspaces -E-n, E-n+2, ■ ■ ■ , En^i form a single orbit under T^. 
Likewise, using o-^crf yields another orbit under formed by the subspaces En-2, ■ ■ ■ , E-n+i- These 
two T^-orbits are clearly related by a global flip Vi and thus form a single orbit under T. 
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This completes the classification of the T-orbits: they coincide with the sectors in (j3.9p in the open 
cases and in the even periodic case, whereas in the odd periodic case, the full space forms a single T-orbit. 
Consequently, the T^-orbits are given by the T^-sectors in (j3.10p in the open and even periodic cases, 
because preserves the value of v in these cases. In the odd periodic case, the above arguments show 
that the two T^-orbits are given by 

n n — 1 

£+= E,, E,. (3.16) 

t)=— ?i+l,by2 v=—n,hy2 

In what follows, we determine the conformal content of the various sectors discussed above. 



IV. CONFORMAL SPECTRA AND PARTITION FUNCTIONS 
A. Diagonalisation of the transfer matrix 

For the two boundary conditions, periodic or open, it is easy to see that T is symmetric and therefore 
diagonalisable. For a € M, T is furthermore real. The explicit diagonalisation may be achieved through a 
Jordan- Wigner transformation. Because the non-periodic case is slightly simpler, and since the periodic 
case was treated by Lieb in [l|, we will briefly outline the diagonalisation procedure in the open case. In 
view of the results of the previous section, it is natural to diagonalise (diagonalising T itself is possible 
by the same technique, but involves unnecessary complications), 

7V-1 N-1 

= V^VI = exp (a ^ c^r^~hi) (« ^t^t+i)- (oP^n case) (4.1) 

2=1 i=l 

We start by finding the spectrum of T^. From this, in Section |Vl we form the eigenbasis for V. 
Following Lieb, we define the following fermionic operators, 

i-1 i-1 

Q = (-1)^"' (H ^fc) , Cl = (-l)-i (H of, {l<i< N). (4.2) 

k=l k=l 

The algebra satisfied by these operators follows easily from the algebra of the Pauli matrices, 

{C„ C,} = {Cl C]} = 0, {Q, C]} = U,^, . (4.3) 
From (j4.2p . we readily obtain Ci Cj+i = —a~ cr~^-^^, and thus 

N-l N-1 

= exp {-a ^ aCi+i) exp (a Cjcl^). (4.4) 

i=l i=l 

The quadratic forms in the exponentials can be block-diagonalised by the following Fourier transforma- 
tion. 



, . I X N 



k=l k=l 
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with e^iT = (-i)^^ for convenience. The transformation being unitary, the fermions Vk^Vk satisfy the same 
algebra as the fermionic operators C-,Cj. Under this transformation, acquires the factorized form 



= (§) Ak, ^fc = exp|2acosgfc %??;v+i-fc} exp|2acosgfc ryjY_^^_^?7j|, gfc = ■ (4.6) 

A:=l 

We note that for A'' odd and k = the block Ak = I2 is the identity (cos qk = 0) and acts in the two- 

dimensional space generated by |0) and J^lv+i |0). For k < y, each block Ak acts in the four-dimensional 

2 

Fock space generated by |0), ?/].|0), f/jy+i-fcl^) ^^"^ ^I^]v+i-fcl^)- -'-'^ cases, the full space has dimension 
2^. 

The diagonalisation of the blocks Ak is straightforward and yields the following eigenvalues. 



Afc = 1 , 1 , + a2 cos2 Qk + a cos qkf , [\/l + cos^ - a cos qk]'^ , (A; < ^) (4-7) 
Afc = 1 , 1 . {k = ^en) (4.8) 



In (j4.7p . the two eigenvalues 1 correspond to the odd (fermonic) eigenvectors ??|.|0) and ^?]v+i-fel^)' 
while the other two eigenvectors are suitable linear combinations of the even (bosonic) states |0) and 

vlvl+i^kl^)- 

It follows that the eigenvalues of are of the form 



LfJ 



A = Afe (4.9) 



k=l 

with the possible values of each Afc as listed above. For A'' odd, each such eigenvalue must be counted 

twice due to the presence of the block A n+i = l2- For both parities, the eigenvalues can be conveniently 

2 

written as 

N-l 

-1 / 1 J- rv'^ Gin-^ Til -L Gin ni tit — 

2(A^ + 1) ' 



nr r. • i2(i-efc-/ifc) -rrk 

yi + a^sm^Pfc + asmpfc , pk = r,/^ , -.n > (4-10) 



fe = JV— 1 mod 2 



where for each k, the numbers Sk and Hk independently take the values or 1. This convention differs 
from the one used in [11] on critical dense polymers. 

We finish this section by computing the conformal limit of this spectrum. The energy eigenvalues are 
obtained from the above eigenvalues A as 

E = -^logX. (4.11) 
We consider the cases A'^ even and N odd separately. 
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B. Conformal spectra for N even 

Here we consider the case even. Since sinp^ > for all k, the (non-degenerate) maximal eigenvalue 
Amax is obtained by choosing all Sk = fik = 0- It follows that the ground-state energy is given by 

^ N-l N-l 

Eo = -- log Amax = - ^ log Y 1 + a2 sin^ Pk + a smpk^ = - ^ arcsinh(a sin ^{n+i) '^ 

fc=l,odd fc=l,odd 

N + 1 ^ , , , 1 , avr ,^ , 

= / dt arcsmhiasmn H — arcsmha h • ■ ■ , (4.12) 

TT Jo ^ ' 2 24N ' ^ ^ 

where the last expression indicates the first few terms in the asymptotic expansion of Eq with respect to 
the system size N. Using the convention 

^0 = A^/bulk + /bdy + O(^) (4.13) 

to identify the free energies, it follows that the bulk free energy per site (half a dimer) is 

. I , . , . . ^ 1 1 arctan(t) i , , ^ ,^ 

/bulk = / dt arcsmh(Q! sm t j = / at = — X2(i«j> (4.14) 

see Appendix [BJ For a = 1, the Legendre chi function in (j4.14p reduces to Catalan's constant, 

- ix2i}) = G = 0.915965594177 .... (4.15) 

Likewise, the (total) boundary free energy is given by 

1 1 /"^ 1 i 

/bdy = - arcsinh a / df arcsinh(asint) = - arcsinha H X2^^^)- (4-16) 

2 vr Jo 2 vr 



Interpreting the -^-correction in (I4.12|) as 

+ (4.17) 



N 



24iV N \ 24 
we read off the effective central charge 

Ceff = c-24Ar" = l (4.18) 

where Ag™"^ is the conformal weight of the ground state for even. 
The asymptotic expansion for large of the excited levels reads 

N-l 

E-Eo= Yl + IJ-k) arcsinh(a sin 2^^) = ^ Yl {i^k + IJ-k) k + . . . . (4.19) 

fc=l,odd A:>l,odd 

The excited energy levels are thus of the form (by redefining e^, //^ — )• ej, fij for k = 2j — 1) 

Ep = Eo + ^p, p=Y i^j + Mi) U - I)- (4.20) 
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The generating function for the level degeneracies dp is then given by 



(4.21) 



where r/(z) is the Dedekind eta function ()A.7p . Here and in the following, 9j is the j-th standard Jacobi 
theta function given in Appendix [Xj We obtain the conformal spectrum generating function, 



Ziq) = E e~ 



EM , „HZlM„ (^3{q) 

= e 24]v N dpC N P — 

p=0 



e N , 



(4.22) 



Up to a renormalisation, it is equal to Tr (T^) 2 and may be interpreted as the universal partition 
function of the dimer model on a cylinder of even height and even perimeter M. The result ()4.22p is 



well known {371 . 
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C. Conformal spectra for N odd 

The case A^ odd is very similar to the even case. The ground-state energy, now doubly degenerate, is 

N-l 

Eq = — arcsinh(a sin 2(n+i) ) 

k=2, even 



7r 

f 2 1 OiTT 

/ dt arcsinh(a sin t) -\ — arcsinh a H — : + . . . . (4.23) 

JO 2 12A' 



_ N + 1 f2 . _ ^ . . 1 . , .an 

vr 

The bulk and boundary free energies are as for A^ even, (|4.14p and (|4.16p . while the ground-state conformal 



weight for A^ odd is related to the one for A^ even by 

Af'^ = Ag™^ + -. (4.24) 
8 

The values of the excited levels are given by 

N-l 

E - Eq= ^ {ek + Hk) arcsinh(a sin ^^n+i) ) ^ 2N ^ i^k + l^k) k + . . . , (4.25) 

fc=2,even fc>2,even 

each such eigenvalue being counted twice. The excited energy levels are thus of the form (by redefining 

i^k,fJ'k £j,fJ'j for k = 2j) 

Ep = Eo + ^p, p = Y^ (e, + (4.26) 
It follows that the generating function for the degeneracies is 

5^ dp = 2 n (1 + 2g^' + q'n = 2 n (1 + Q' f = 4y • (4.27) 

p=0 j=l j=l ^^'i' 

The spectrum generating function, identified with the universal partition function on a cylinder of odd 
height A^ and even perimeter M, becomes 



37 



Z(g) = ^ e"^^^ = e-w2_] dpe"— ^ = q = e- — . (4.28) 



p=0 
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V. SECTORS AND REFINED SPECTRA 

Using the definition (|4.2p of tlie auxiliary fermions Cj,Cj in terms of Pauli matrices, and then the 
Fourier transform (j4.5p to the fermions ?7a:,?7|., the operator V can be written as 

1 / i\Af ^ -] ( -\\N ^ 

V = '-^^ + E (-1)^' C]C, = + i-- 5: (-l)'^ vi ri.^.-,- (5.1) 

j=i k=i 

Now, recalhng that factorizes as 

= ^fc, ^fc = exp |2a cosgfc %?/Ar+i-A:} exp|2a cosgfc r/]y^-^_^7?||, (5.2) 

k=l 

we see that the terms in the sum (|5.1|) should be grouped in pairs to act on the same spaces as the 
blocks ^fc. This allows us to examine each block separately. As seen in the following, the action depends 
slightly on the parity of N. Before examining the spectrum in the various sectors, we discuss the relation 
between the variation index and the fermionic parity. 

A. Variation index and the fermion number 

The Jordan- Wigner transformation converts the arrow configuration space into a fermionic Fock space 
whose states are built by applying the discrete fermions r/^ on the vacuum state |0). The vacuum, being 
annihilated by all r/^, hence by all Cj, is proportional to ••• |). The full space may be split into 
an even sector and an odd sector, made up respectively of states containing an even or odd number of 
fermionic excitations. From the explicit expression (12. 6p of the transfer matrix, and the relations (j4.2p 
between the Pauli matrices and the fermions, we see that T leaves the even and odd sectors invariant if 

is even, but interchanges them if is odd. The same is true in the periodic case [li]. It follows that 
the even and odd sectors are left invariant by T^. The above expression ()5.2p makes this manifest in the 
open case. 

In all cases, open or periodic, the fermion number operator, given by 

N N N 

k=l i=l i=l 

is closely related to the variation index operator V. In particular, the fermionic parity (—1)^ is related 
to the variation index by 

f (-1)^ if N is even, 
[ (-1)^-2 if AT is odd. 

In the periodic case with N odd, the two invariant sectors iS+ and £- coincide with the even and 
odd subspaces, although the precise identification depends on N. Since the vacuum state, which is an 
even state, has variation index v = we obtain from (|3.16p that the even subspace is equal to for 

= 1 mod 4, and to £- for A^ = 3 mod 4. 
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B. Refined spectra for N even 

For even, the operator V is rewriten as 

E. E. 

V = ^ (-1)^ [4 VN+i-k - vUl-k Vk] = E (5-5) 

k=l k=l 

thereby mtroducing the operators Vk as the indicated summands. For every k, the operator Vk is di- 
agonahsable and commutes with the corresponding block Af.. It follows that the two non-degenerate 
eigenvectors of Ak are eigenvectors of Vk as well. Indeed, we find 

Vk\0)=VkvlvUi-k\0)=0. (5.6) 

Also, it is recalled that the two eigenvectors of A^ corresponding to the degenerate eigenvalue 1 are r/]. |0) 
and VN+i^k -'■^ basis of these two vectors, Vk is represented by the matrix 

Vfc=i^-i(-l)'^ (^^^ , (5.7) 

with eigenvalues Vk = ±i^(— l)'^ G {—1,1}. The degenerate Ak eigenspace of eigenvalue 1 is therefore 
spanned by two Vk eigenvectors with the distinct eigenvalues +1 and —1. 
Returning to the spectrum of Ak given earlier, namely 



2(l-£fc-^Jfc) 

1 + a'^ sin pk + asinpk , (5-8) 



we see that the two states with + /i^ = 0, 2 have Vk = whereas those with e/j + /i^ = 1 have Vk = ±1. 
Introducing an extra parameter y to keep track of the eigenvalues of V, we obtain the following generating 
function for the degeneracies dp^y at level p and eigenvalue u of V, 

p=o vez j=i 

Taking the ground-state energy into account, the corresponding conformal partition function reads 

Z{q- y) = Tr (T^^y^) =q-^^Y.Y. ^P,^ = <1 = ■ (5-10) 

p=0 dSZ 

By expanding ^3 in a power series in y and comparing with 

Z{q;y) = Y,ZM)y'', (5-11) 

II 

we obtain the partition functions pertaining to the individual T^-sectors, 

it 

Z.(g) = ^, v€Z. (5.12) 

It is recalled that this derivation of the sector partition functions (|5.12p is based on T^, so M is assumed 
even. 
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C. Refined spectra for N odd 

For odd, the operator V in ()5.ip may be written as 

N~l 
2 

+ i^-i J2i-'^)'\vUN+i-k + vWi^k^k ^Vm^ + Y, ^k, (5.13) 



V 



n t 

- — Vn+i V n+1 

^ ~2~ 2 



iV-1 iV-1 
2 2 

J 

k=l " " A:=l 



thereby introducing the (diagonalisable) operators Vk- For every k < the operator Vk acts in the 

same space as the block with which it commutes. Both A n+i = I2 and V n+i are two-dimensional. 

2 2 

For A; < ^^y^, as before, we find that the non-degenerate eigenvectors of have Vk = 0. Moreover, 
on the two degenerate states rjl |0) and ??]v+i_fc |0), Vk is represented by 

i^-'i-l)'^ (5.14) 

with eigenvalues Vk = =bi^~-^(— l)'^ G 1}- We reach the same conclusion as in the N even case: the 
two eigenvectors of Ak with eigenvalue 1 have Vk = ±1. 

In the remaining two-dimensional block, the operator Vjv+i = ^ — r/lr+i Vn+i takes respectively the 

2 ~2~ 

values ^ and — ^ on the states |0) and ?7jv+i|0), themselves eigenvectors of A n+i with eigenvalue 1. 

While in the previous section, this doubly degenerate eigenvalue 1 of A n+i caused the duplication of all 

2 

eigenvalues of T^, it now gives an extra contribution vn+i = i to half the spectrum, and v n+i = — i to 

2 ^ 2 ^ 

the other half. 

The corresponding generating function for the degeneracies dp^^ reads 

f; E dp,,qPy- = iy^+y-^)f[{l + {y + y-')q^ + q^^)=q-^2^-2^, (5.15) 



and leads to the following conformal partition function. 



Z(g;,)=Tr(r*V)=<Z^E ^ dp,l'y'="-^^ g = e'^. (5.16) 

Prom the power expansion of 62 with respect to y, we obtain the partition functions pertaining to the 
individual T^-sectors, 



^.(9) = ^, veZ + l, (5.17) 



again assuming M even. 

Strikingly, for every N , even or odd, we obtain the unified conclusion that in the conformal limit, the 
sectors are labelled by (half-)integers v G Z + and the partition function of a single sector is equal to 
the character of a single Verma module, 

(7 2 24 

Z,{q) = Z^,{q) = . (5.18) 

n (1 - ^'") 
m=l 
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However, as discussed in the following, we do not claim that the underlying representations necessarily 
correspond to Verma modules, only that the characters are as given in (jS.lSp . 

VI. CONFORMAL FIELD THEORY 

It is natural to ask what conformal representations arise in the continuum scaling limit of the various 
sectors. As in other comparable studies, the answer is not immediate from the knowledge of the transfer 
matrix alone. In particular, the possibility of non-trivial Jordan blocks in the conformal limit is a 
subtle issue. The transfer matrix under investigation here is diagonalisable. However, many eigenvalues, 
which are distinct for finite system sizes, become degenerate in the conformal spectrum, thereby making 
room for non-trivial Jordan blocks. Thus, the diagonalisability of a finite-size transfer matrix does not 
necessarily preclude the existence of such Jordan blocks in the conformal representations, although no 
example of this kind seems to be known^. We intend to address this possibility for the dimer model 
elsewhere. In addition, as is well-known, the value of the central charge c is crucial to characterise the 
associated conformal field theory and the structure of the various representations, but this value cannot 
be determined unequivocally from the spectrum alone. 

Here we will make the ansatz that the dimer model has 

c=-2. (6.1) 

This value is borne out by a number of independent observations and results: (i) the dependence of the 
effective central charge on the parity of the width of a strip is easily and naturally explained from that 
particular value of c [6]; (ii) there is a well-known equivalence between the dimer model in a rectangular 
geometry, the uniform spanning tree model and the Abelian sandpile model, the latter providing over- 



whelming evidence that c 



27 



311]: (iii) the partition function on a rectangle with odd side lengths 
and one corner site removed is equal to ri{q) in agreement with the general result [ri{q)\~'^^'^ , see 42]; and 



(iv) the finite-size corrections to the conformal spectrum can be reproduced in perturbation theory for 
c = —2, but appear unreproducible for c = 1 (as we will discuss elsewhere). Moreover, we will see that 
the conformal spectrum of the dimer model is very closely related to that in the critical dense polymer 
model, also believed to have c = —2, see [ll, [3], for example. 

Notwithstanding, based on different observations, there are also strong arguments favouring the other 
widely proposed value, namely c = 1: (i) the so-called height function, defined on dimer configurations. 



converges in the scaling limit to a Gaussian free field 



431 ]: and (ii) in the generalised model allowing 
le monomer correlators are those of free 



for monomers in addition to dimers, it has been shown that t 

complex fermions, both on a boundary {2^ and in the bulk 21], following earlier observations 17l4l9l]. 



In the critical dense polymer model, for instance, the non-trivial Jordan blocks are present for finite system sizes [Ill,l40l.l41| 
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Taken together, these results and observations seem to indicate that the pure dimer model has c = —2, 
while the monomer-dimer model has c = 1. In this scenario, the height function should be viewed as 
belonging to the monomer-dimer model. 

Now, an immediate consequence of assuming c = —2 is that the parity-dependent ground-state con- 
formal weights are 



A even 1 a odd ri 



■2) 



while the representation associated with the partition function (q) in the sector v has conformal weight 



A(v) 



4t>^ - 1 



V G 



3) 



8 ' ' ~ ^ 

These weights are integer if N is odd (for which v G Z + i), but equal to — | modulo an integer or a 
half-integer if is even (for which f G Z). They can be expressed by the Kac formula for c = —2 as 

(2r-s)2-l 



^r,s 



r, s G N, 



(6.4) 



where r, s are Kac labels. The corresponding extended Kac table is shown in Figure [2j Specifically, the 
set of conformal weights in (j6.3p are covered exactly once by setting 

, , s 



\v\=r - -, r G N, s G {1,2}. 



(6.5) 



This corresponds to the two lower rows in the Kac table. Following [1^, [32], and without claiming any 
super conformal symmetry in the dimer model, we refer to them as the Neveu-Schwarz (s = 1) and 
Ramond {s = 2) sectors. Alternatively, the conformal weights in (j6.3p are covered exactly once by setting 



s 

r, 

2 



r = l, s = 2,3 



(6.6) 



noting that Ai i = Ai 3 = 0. As discussed below, this labelling offers a direct link to the characterisation 



Hi 



of a particular class of boundary conditions in the critical dense polymer model 

With the ansatz (j6.ip . none of the representations associated with the f -sectors is irreducible. Indeed, 
in the case of open boundary conditions, we have the decompositions 



rj{q) 



rj{q) 



ch 4(].|+2,)2-i (g) = ch^(|t,|+2j)(9)' 

j=0 * j=0 

00 00 

Y ch 4(|,|+,)2_i (g) = Y ch^(|H+i)(9), 



V G 



2' 



(Ramond) (6.7) 
(Neveu — Schwarz) (6.8) 



j=o j=o 
in terms of irreducible c = — 2 characters. 



chA(m)(9) 



g^(l_g2H+2) 

viQ) 

g4(l-gl™l+l) 



m G 



m G Z+ i. 



(6.9) 
(6.10) 
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FIG. 2: Extended Kac table for c : 
sectors. 



-2. The two lower rows label the Neveu-Schwarz (s — 1) and Ramond (s — 2) 



also use the alternative notation introduced in 



Here we have adopted the notation ch.^{q) used in 11[ to indicate irreducible c = —2 characters. We will 



ch (q) = ch^^^(g), 



(6.11) 



utilising the Kac labels r, s. 

The states in the representations underlying the partition functions in (|6.7p and (j6.8p correspond 



to eigenvectors of the transfer matrix. The combinatorics encapsulating this correspondence is closely 
related to the physical combinatorics of the critical dense polymer model as described in [ll|, [l^. This 
intriguing observation is outlined below and depends on the parity of N. 

Here we note that the partition functions in (j6.7p and ()6.8p are characters of Verma modules of 
highest weight A{v). We stress that this does not necessarily imply that the underlying modules are 
Verma modules, only that they have Verma module characters. In fact, it is beyond the scope of the 
present work to characterise the specific module structures. 

VII. OPEN BOUNDARY CONDITIONS 
A. Physical combinatorics in the Ramond sector: TV even 



From Section HYI with c = — 2 and Aq™" = —k^ the conformal energy levels for N even read 



^ ^ an X - 1 avr r c 1 

E = Eo-\ > (]-h)6i = — <^ 

" V ^ 2^ J ^ \ 24 8 

i>i i>i 



(7.1) 



21 



where 6j = ej + fij = 0, 1, 1 or 2. The results of Section |V] show that 6j = 0, 2 do not contribute to the 
value of the variation index v, while 6j = 1 contributes +1 or —1 to v. Our convention will be that Sj = 1 
contributes —1 to v, while ^ij = 1 contributes +1, so that the value of of a state specified by the set 
{eii/^i}i>i is given by 

i>i 

For finite excitations, the upper bound in this summation is given by n = y- 



From (j6.7p . the partition functions of the first few T;-sectors read 

Zo{q) =ch_i(q) + chi5(g) + ch63(Q)+chi43(g) + ..., (7.3) 

8 8 8 8 

= ch3(g)+ ch35(g)+ch 99(g) +ch 195(g) + ..., Z_i(g) = Zi(g), (7.4) 

8 8 8 8 

^2(9) = chi5 (q) + ch63 (g) + chi43 (g) + ch255 (g) + . . . , Z_2{q) = Z2(q). (7.5) 

8 8 8 8 

The reason for indicating some of the characters in boldface is described in the following. 



There is a striking similarity 
polymer (CDP) model studied in 



between the spectrum (j7.ip and the spectrum of the critical dense 
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12l |. In fact, the conformal energies (j7.ip are exactly those of the 
CDP model on a strip of odd width (rather than even as here) provided the values {ej, fij}j>i obey certain 
selection rules. In other words, the CDP and dimer conformal spectra differ only by the degeneracies 
of the energy levels. This observation also applies to the Neveu-Schwarz sectors where v G Z + ^. In 
addition, a significant part of the CDP conformal spectrum can also be partitioned into sectors^, labelled 
by a non-negative integer £ counting the number of so-called defects. This defect number is related to 
the dimer sectors with v > through 

i = 2v + l, vG^No-^. (7.6) 

This relation is intimately associated with ()6.6p . The equations (I7.3p - (l7.5p show the conformal content of 
the first few Ramond sectors in the dimer model; the characters in boldface are those which appear in the 
CDP model for the corresponding defect number. For instance, the irreducible characters of conformal 
weights — |, |, ^ appear respectively once, twice and three times in the dimer model, but only once in 
the CDP model as i varies over the odd positive integers. 

The question still remains of how to associate the eigenvalues in (j7.ip to the various representations. 
This may be answered in a way reminiscent of the corresponding results in the CDP model, but in a 



^ Although we use the word "sector" in the context of the CDP model, as in [ij], it does not have the same meaning as 
in the dimer model. As discussed above, the variation index « is a good quantum number in the dimer model. In the 
CDP model, on the other hand, the defect number £ is not a good quantum number since the corresponding subspaces 
are not left invariant by the finite-size transfer matrices. Instead, the CDP transfer matrices are upper block triangular 
with respect to £ [ill ). 
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i = 7 
j = e 

j = 5 
j = 4 

i = 3 
i = 2 
i = i 



FIG. 3: Two-column configuration characterised by non- vanishing occupation numbers £21 £31 £67 /^ii /^2j fJ■4^ fJ-s- This 
is an example of a non-admissible configuration. 



somewhat surprising manner as it is based on the CDP model defined not only on the strip 11[ , but also 
on the cylinder 

Following [n I , we thus associate an eigenvalue of the form ()7.ip with a two-column (double-column) 
semi-infinite configuration as in Figure [3l where the levels, labelled by j > 1 from below, contain two 
sites, empty (white) or occupied (grey). For a given set {ej;/^j}j>i; the left site at level j is occupied 
if Ej = 1 but empty if ej = 0, while the right site at level j is occupied if fij = 1 but empty if fij = 0. 
Finite excitations of the form (j7.ip correspond to configurations with only finitely many occupants. As 
indicated, we occasionally refer to an occupied site as an occupant. 

For finite system size N = 2n and variation index v, a natural finitization of the partition function 
Zy(q) is now given by 



c I A even 



E 

1"— two 
columns 



(7.7) 



where the sum is over all two-column configurations of maximal height n and variation index v, and 
where the elementary excitation energy of an occupant at level j follows from (|7.ip and is given by 



Ej=j-l 



(7.8) 



The two-column configuration in Figure O for example, has energy excitation E — Eq = ^ ■ ^ and 
belongs to the sector v = 1. The non-zero values of 6j = Ej + fij are Si = 63 = 6^ = 65 = 6g = 1 and 
62 = 2 implying that this state is degenerate with 31 other states with the same energy and values of 5j. 
The corresponding two-column configurations are obtained from the present one by moving occupants 
from left to right, or from right to left, at the levels containing a single occupant. These 31 other states 
do not all belong to the same sector, though, since moving an occupant horizontally changes the value of 
V. Each of these states will have a variation index v between —5 and -|-5, since there are 5 distinct levels 
with a single occupied site. This also accounts for the total of 2^ = 32 related configurations above. 
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According to [11], a two-column configuration with a finite number of occupants (corresponding to a 
finite energy excitation) is admissible or non-admissible^. An admissible configuration is such that if one 
draws a horizontal line at any height, the number of occupants on or above the line in the left column 
must be smaller than or equal to the number of occupants on or above the line in the right column. 
Equivalently, admissible configurations have values {ej,iJ,j} such that 

ifJ-j - Sj) > 0, Vk G N. (7.9) 

In particular, they all have a smaller or equal number of occupants in the left column than in the right 
one, and consequently have values v > 0. It is noted, though, that a two-column with v >0 need not be 
admissible, as illustrated in Figure [3l 

The set of admissible configurations of maximal height h (i.e. no occupants above level h) and given 
variation index v may be decomposed into a disjoint union of subsets, 

h—v 
m=0 

where contains the admissible configurations of maximal height h with exactly m occupants in 

the left column and exactly m + v occupants in the right column. Thus, Ai^^ is the set of admissible 
configurations which have an excess of v occupants in the right column over the left column. In the 
continuum scaling limit, where there is no upper bound on h, the set of such admissible configurations 
becomes 

X = lim Ai^\ (7.11) 



For the CDP model, Pearce and Rasmussen 



111 ] conjectured a set of selection rules positing that the 



CDP conformal spectrum on a strip of odd width is given by (j7.ip for all those values {ej, which define 



44]. 



admissible two-column configurations. This conjecture has been recently proven by Morin-Duchesne 
In the continuum scaling limit, the combinatorics of admissible configurations is such that the set Av 
corresponds to conformal states exactly matching the degeneracies of a single irreducible representation 
of highest weight A(v). That is, the CDP partition function in the sector with defect number i = 2v + 1 
is given by 

ZCDP(g) = Zg1\{q) = ch^iq), V G Nq (7.12) 

8 

corresponding to a boldfaced term in the expansions given above in (j7.3p - ()7.5p . In other words, the 
boldfaced terms comprise the contributions from the admissible configurations. 



We have chosen to adopt this terminology from the work on the CDP model where it is very appropriate. In the dimer 
model, on the other hand, non-admissible configurations actually carry physical information about states in sub-leading 
characters. A more befitting name for these configurations could therefore be subordinate configurations. 
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FIG. 4: Ramond sectors (N even): Two-column configurations with minimal energy excitation (I7.14p for fixed 
quantum number v given by the excess of blue (right) over red (left). 



In the dimer model for A'" even, as already indicated in (j7.7p . all two-column configurations of 

maximal height : n = ^ (Ramond) (7-13) 

must be included, even though many of them in general are non-admissible. They nevertheless carry 
physical information and are accounted for by a combinatorial enumeration reminiscent of the physical 
combinatorics of the CDP model on the cylinder [12], likewise for N even, as outlined in the following. 

The two-column configurations with variation index v are generated combinatorially by starting with 
the minimum-energy configuration for given v as shown in Figure HI The energy excitation of such a 
ground state is 



\v\ 

E{v) =Y,{j-\) = \v\ V E Z, 

3=1 

while the corresponding conformal weight is 

A(^;) = Ar" + i^(t') = A|„|+i,2 . 



(7.14) 



(7.15) 



The generating function for the full excitation spectrum for finite system size N is given by the finitized 
version of (15.91). 



n E 

i=l £j,AtjG{0,l} 



N-l 



A:=0 



N 
n — V 



y , 



(7.16) 



from which it follows that the finite-size generating function for fixed v is incorporated in the g-binomial 



-Eiv) 



E 

ii-two 
columns 



n — V 



(7.17) 



Indeed, the elementary excitations (of energy 1) are generated by either inserting a left-right pair of 
occupants at level j = 1 or promoting an occupant at level j to level j + 1 in the same column. These are 
the only operations increasing the energy by 1 while preserving v. The generation of the corresponding 
spectrum and its link to (j7.17p are illustrated in Figure [5l We stress that, even though 

(7.18) 



N 
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n — V 


q 


n + V 
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FIG. 5: Ramond sectors {N even): Combinatorial enumeration by two-column configurations of the g-binomial 
[n-ti] "[2] ~ '^q^i^^^^ ■ The excess of blue (right) over red (left) is given by the quantum number v = 1. 
The lowest energy configuration has energy E{v) — = ^- Of the 15 configurations appearing here, only the 
top one is non- admissible. 

as q-polynomials, they have different combinatorial interpretations because they have different quantum 
numbers: v and — f , respectively. The combinatorial interpretations obviously agree in the exceptional 
case where v = Q. 

This combinatorial enumeration immediately leads to the finitized partition function 



"1+1.2 



N 
n — V 



\v\ < n 



N 



(7.19) 



and our next objective is to write this in terms of finitized irreducible c = — 2 characters [111 . 11/ 

In the Ramond sector, the relevant such characters are the ones associated with the second row of the 
extended Kac table in Figure [21 and for which the finitization is with respect to an odd strip in the CDP 
model. Following {l^ . we thus have 



(N+l) 
r,2 



where the finitized character ch^^^"^^ {q) 



n or n+l 
r=|t;| + l,by2 

di\"^^^\q) for s = 2 is given by 



(7.20) 



ch, 



(N+l) 

T,2 



(9) 



,2r 



,Ar+2 



N + 2 
n — r + 1 



(7.21) 



and where the upper summation limit depends on the parity oi n — v. Since ch^^^^\q) vanishes for 
r > n + 1, the specification of this upper limit is redundant, as it could be replaced by 00. In the 
continuum scaling limit — t- 00, we readily recover the decomposition (j6.7p . 



limZi^)(g)= V cK2{q) = Zy{q). (Ramond) (7.22) 

Af— i>oo ^ — ' ' 
r=|i;| + l,by2 

It is recalled that only the first term c]i^v\+i,2{q) (|6.7p and (|7.22p is related to the CDP model 
through ()7.6p . as indicated in (j7.3p - ()7.5p . Likewise for finite system sizes, only some eigenvalues con- 
tribute to ch|^^j^2(Q)) namely the ones corresponding to admissible two-column configurations. Here we 
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V ^ L Wadm = O 

FIG. 6: Illustration of how our prescription turns a non-admissible two-column configuration into an admissible 
one. The non- admissibility measure of the initial configuration is t = I, and its single non- admissibility label is 
Ki — 3. The variation index is increased by a multiple of 2, here from v — 1 to fadm = w + 2< = 3. This is the 
maximal attainable value, cf. (I7.3ip . The resulting configuration is the ground-state configuration for Wadm — 3. In 
accordance with (|7.25p . its energy excitation satisfies i?(fadm) = E{v) + 4, where E{v) ~ E{1) ~ i. 

present a combinatorial prescription which, for given v and n, associates every corresponding two-column 
configuration with a state in one of the finitized characters appearing in the decomposition (j7.20p . giving 
this decomposition a clear combinatorial content. 

Before giving the general prescription, let us examine (q) featured in Figure El From ()7.19p . we 
have 



whereas (|7.20p yields 

^ {q) = ch2(g) + ch2(g) = (1 + g + 2g2 + 2q^ + 2g^ + 2q^ + 2q^ + + q^) + ^-^+^^,2 _ (7,24) 
Consistency requires 

A2,2 + 4 = A4,2, (7.25) 

which is easily verified. The only two-column configuration in Figure [5] not contributing to ch2 2(g) is the 

single non- admissible configuration appearing at the top. The announced combinatorial prescription turns 

(7) 

it into an admissible configuration with v = 3 contributing to ch4 2(9) by moving the single occupant in 
the left column to the right column, as illustrated in Figure [6j This way of accounting for the sub- leading 
finitized characters in (j7.20p extends to general Z^^\q) as described in the following. 

First, as a measure of 'non-admissibility' of a finite two-column configuration, we introduce 



r-l-A2 



(7.23) 



t 



(7.26) 



According to (j7.9p . an admissible configuration is thus characterised by a vanishing non-admissibility 
measure, t = 0. 

Let C/^ denote the set of two-column configurations of maximal height h and variation index v. The 
cardinality of this set is 

h — V , 



(7.27) 
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FIG. 7: Prescription for turning a non- admissible two-column configuration into an admissible one. The initial 
configuration has non-admissibility labels ki, H2, and variation index v, while the resulting admissible 

configuration has variation index fadm = v + 2t. 



For a configuration in J7^, the non-admissibility measure t satisfies 



max{0, -v} < t < traa.x, t„ 



h — V 



(7.28) 



Now, consider a two-column configuration in with non-admissibility measure t. Associated to this 
configuration and for every r = 1,2, ... ,i, we introduce the non-admissibility label as the maximum 



value of /t G {1, 2, . . . , /i} for which J2j>Kif^j ~ ^j) ~ That 



IS, 



Kt = max 



{k G {1, 2, . . . , h}; - ej) = -r}- 



(7.29) 



By construction, such a configuration must be of the form indicated to the left in Figure [71 The t 
occupants at levels k^, t = 1,2, ... ,t, all appear in the left column with unoccupied sites to their right. 
In addition, the parts of the configuration above the level ki and in-between the levels k^- and Kr+i for 
every r = 1, 2, . . . , t — 1, make up admissible sub-configurations with zero variation index; the part of the 
configuration comprising the lowest — 1 levels makes up an admissible subconfiguration with variation 
index v + t. 

Our prescription to associate an admissible configuration to a non-admissible configuration of the 
previous type is simply to move the t occupants at levels ki, K2, ■ ■ ■ , Kt from the left to the right column 
(as in Figure [7] and illustrated in Figure [6|, resulting in a configuration with variation index 



t'adm =V-\-2t. 



(7.30) 



It remains to be demonstrated that the admissible configurations obtained in this way exactly corre- 
spond to the ones associated with the finitized characters appearing in the decomposition ()7.20p . We do 
this by establishing that the prescription yields a bijective map 



n— 1 orn 

^ U <L 

V!idui=\v\,hy2 



(7.31) 
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where the upper hmit depends on the parity of n — v. The bijectivity of this prescription map follows by 
observing (i) that the cardinality of the union in (j7.3ip equals the one in ()7.27p for h = n; and (ii) that 

(n) 

every admissible two-column configuration in ^4^,^^^^ with variation index written as Wadm = v + 2t, where 
t,v + t > 0, is of the form indicated to the right in Figure [71 

The range of the (disjoint) union index in (j7.3ip is dictated by the possible values for fadm following 
from (j7.28p and (j7.30p . These values are readily seen to be the ones labelling the summation in the 
decomposition (|7.20p . reflecting that the admissible configurations in ^4^^]^^ are the ones contributing to 
the finitized character ch^^"*^^]^ 2(9) (|7.20p . Each of the finitized characters in (j7.20p therefore appears 
as the partition function for a fixed value of v and a fixed non- admissibility measure t > 0, 

C(5) = ch[;i^k2(«)- (7-32) 



It is noted that the generalised Catalan number 



<:;i.M(i)=u-M-2,j-i„-i«i-a-2j- <^-^'> 

defined by evaluating the corresponding finitized character at q = 1, gives the number of two-column 
configurations with variation index v and non-admissibility measure t. 

The finitized full partition function is obtained by ignoring the separation into u-sectors, and is given 

by 

n n+1 



N \ / N 



Z^''\cD= E Zf\q) = Y.rol.Z^'\q). (7.34) 

r=l 



This shows that the finitized character chj,^"*"^^ (g) appears exactly r times in the corresponding full 
partition function provided r < n + 1, and thus explains the noted appearance (once, twice and three 
times) of the irreducible characters of conformal weights Ai^2 = — |, ^2,2 = | and A2,3 = ^ in (j7.3p - (j7.5p . 

We conclude this analysis of the Ramond sector by reconsidering the full partition function in the 
continuum scaling limit, that is, 



Z{q) = lim zW(g) = Vz.(g) = Vrch,,2(g). (7.35) 

N^oo ^ — ' ^ — ' 

In terms of the W-irreducible characters [3, B| 

X_iiQ) = —rT^o,2iq) = ^ rch^,2(g), Xs (q) = — 7-^^2,2(1) = ^ ^K^iQ), (7-36) 

* ^^^^ re2N-l ^^^^ re2N 

the partition function ()7.35p thus reads 

Z{q)=X-l^{q)+Xl{q). (7.37) 
This is in accordance with ()4.22p since Osi^q) = '&o,2{q) + '&2,2iq)- 
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B. Physical combinatorics in the Neveu-Schwarz sector: TV odd 

It is recalled that, for odd, the variation index v takes one of the + 1 values 

vG[-n,n]n{Z-^), n = f G N - i. (7.38) 
From (I4.25p . with c = — 2 and Aq^'^ = 0, the conformal energy levels for odd follow from 

where 5j = 0, 1, 1 or 2. No contribution to the value of v is made if 6j = 0, 2, while 5j = 1 contributes +1 
or —1 to V. In addition, every energy level is doubly degenerate; one of the paired energy levels gets an 
extra contribution to v equal to the other an extra contribution equal to — ^. For finite excitations, 
we thus introduce the excess parameter 

-I 

so that 

w = v±^, v = w^^, w e {-n + 5, -n + |, . . . , n - i} C Z. (7.41) 

Since < n — ^, we have the unique relation w = v — ^ (resp. w = v + ^)iiv = n (resp. v = —n). For 
even, there was no need to introduce a separate excess parameter since the variation index v itself was 
sufficient to distinguish between degenerate energy levels, and we simply had ()7.2|) . As discussed in the 
following, the excess parameter plays a key role in the physical combinatorics for A^ odd. 
As we did for A^ even, let us write 

5j = Ej + /xj, Ej.iij G {0, 1}, i G N. (7.42) 

To every eigenvalue (j7.39p characterised by {£j>/Uj}j>i, we again associate a two-column configuration 
with occupants in the left or right column at levels j for which Sj = 1 or /Xj = 1, respectively. For finite 
excitations, the maximal height of these two-column configurations is 

maximal height : n — ^ = ^^-^^ (Neveu — Schwarz) (7.43) 

and the excess of occupants in the right column over the left column is given by the excess parameter w. 
The notion of admissibility and the non- admissibility measure and labels are as for A^ even. 
From (|6.8p . the partition functions of the first few T;-sectors read 

= cho(q) + chi(g) + ch3(g) + ch6(g) + chio(g) + + • • • , (7-44) 

2 

Zi (g) = cho(q) + chi(q) + ch3(g) + ch6(g) + chio(g) + chi5(g) + . . . , (7.45) 

2 

Zz_[q) =chi(g) +ch3(q)+ch6(g)+chio(g) + chi5(g)+ch2i(g) + ..., Z^i[q) = Zs{q), (7.46) 

2 2 2 

Zsiq) =chs{q) +ch6(q)+chio(9) + chi5(g)+ch2i(g)+ch28(g) + ..., Z_3{q) = Zr^{q). (7.47) 



30 



i = 3 

J = 2 
J = l 



w 

V 



• o 




O O 




O O 




O O 




o o 




O O 




O • 


• o 




• O 




O O 




O O 




o o 




O • 




O • 


• o 




• o 




• O 




O O 




o • 




O • 




O • 


-3 


-2 


-1 





1 


2 


3 


7 
2 


5 
2 


3 
2 




3 
2 


5 
2 


7 
2 



FIG. 8: Neveu-Schwarz sectors {N odd): Two-column configurations with minimal energy excitation (|7.50p for 
fixed quantum number v. The excess of blue (right) over red (left) is given by the excess parameter w. For these 
minimal energy configurations, w is related to the variation index vhyw = v— ^iiv>0 and w = v + ^ ii v < 0. 

The rationale for writing some of these characters in boldface is as above and is recapitulated below. 

For finite system size N and variation index v, a natural finitization of the partition function Zy{q) is 
given by 



c I A odd 



v-two 
columns 



(7.48) 



where the sum is over all two-column configurations of maximal height n — ^ and excess parameter 
w = V ^ ^, and where the elementary excitation energy of an occupant at level j follows from ()7.39p and 
is given by 



(7.49) 



The two-column configurations with variation index v are generated combinatorially by starting with 
the minimum-energy configuration for given v as shown in Figure El The energy excitation of such a 
ground state is 

E{v) = lv'^-l, ueZ-i, (7.50) 
while the corresponding conformal weight is 

A(^;) = AS^'i + £;(i;) = A,,|_,i^i. (7.51) 
The finitized version of (|5.15p yields the generating function of the full spectrum for finite system size, 



N-l 



i=i £j,fije{o,i} k=o 
and immediately leads to the finitized partition function 



E 



N 
n — V 



N 
n — V 



\v\ < n. 



(7.52) 



(7.53) 



In this case, the elementary excitations (of energy 1) are generated by either inserting a left or right 
occupant at level j = I, provided the excess parameter is w = v-\-^ or w = v—^, respectively, or promoting 



31 



i 


= 3 


O 


o 


j 


= 2 


o 





j 


= 1 


o 


• 



= 1 



o 


o 


o 


• 


o 


o 






o 


o 


• 





• 



o 


• 





o 


o 


o 



2q^ 






• 


o 


o 





• 



o 


o 





• 


• 


• 



2q^ 



o • 




• o 




o • 


• 




• 




• 


o 




• 




• • 








o 




• 




• • 


• • 




• o 




o o 


• 




• 




• 








o • 




o • 




o • 


o 




• 




• • 


• • 




• o 




o o 


3g4 + 


+ 





o 


• 


• 


• 


o 


• 



• 


• 


o 


• 


o 


o 



• 


• 


o 


• 


o 


• 



o • 



2q'^ + 2q' 



,8 



• 


• 




• 


• 


o 


• 




• 


• 


• 


• 




o 


• 



FIG. 9: Neveu-Schwarz sectors {N odd): Combinatorial enumeration by two-column configurations of the q- 
binomial =[2] = S ^i^i _ The excess of blue (right) over red (left) is given by the excess parameter 
u; = 1 or w = 2 related to the quantum number v — n~m — | through w — vzt^. The lowest energy configuration 
has energy E{v) = — i = 1. Of the 21 configurations appearing here, only the top one in the middle column 
is non-admissible. The six configurations appearing at the top of the 3rd, 4th, 5th, 7th, 8th and 9th columns have 
w = 2] the remaining 15 configurations have w = 1. 



an occupant at level j to level j -|- 1 in the same column. These are the only operations increasing the 
energy by 1 while preserving v. The generation of the corresponding spectrum is illustrated in Figure O 
Following [12], we then have 



n+i 



n— i orn-l-i 



AN+l), 



(7.54) 



w=v±\ r=\w\ + l,hy2 



where the finitized characters are given by 

ch(^+l)(<;)=g-^+A.,,. 



1 



n 



N + 1 
r + \ 



The leading character in (I7.54|) appears for r = |v| -|- ^. In the continuum scaling limit 
readily recover the decomposition (j6.8p . 



lim 



{q) = chr,i((?) = Zv{q). (Neveu-Schwarz) 



(7.55) 



00, we 



(7.56) 



r=\v\ + 4 



Admissibility of a two-column configuration requires t = and hence tt; > 0, so with the CDP defect 
number £ > related to the variation index v hj £ = 2v + 1, the characters also present in the CDP 
model, as i varies over the even non-negative integers, are the ones in boldface in (j7.44p - (j7.47p . The 
appearance of two characters in boldface for £ > reflects that the corresponding CDP representation is 
reducible [ll|, 3]. The decomposition of the associated set of admissible configurations of maximal height 



h reads 



h—w 



jih) ^ I I jih) 



(7.57) 



m=0 
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In somewhat sloppy notation, we thus have 

Zi^^) D X"-^) , 4^) D (^1:^^ U ^1;^)) , . G N - i, (7.58) 
whose continuum scahng hmit is indicated by 

z^d (a^_iUA^_^i), ven-l (7.59) 



It is noted that the generahsed Catalan numbers defined by evaluating the finitized characters ()7.55p 
at g = 1, 

verify that 

E ch;.r^(i) = (JJ=4"^(i)- (7.61) 

Separating this sum into the two corresponding w streams as in the double sum in (|7.54|) . we have the 
refined identities 



r=lD-i|+l,by2 Z< r=\v+l\+l,by2 ^ < Z< . 

(7.62) 

Our combinatorial prescription for associating non-admissible configurations with physical states in 
the finitized characters and partition functions extends from the Ramond sector described in Section TVII Al 
to the Neveu-Schwarz sector discussed here. A priori, there appears to be a caveat, though, as moving 
occupants from the left to the right column increases the excess parameter w by a multiple of 2, 

■u^adm = w + 2t, (7.63) 

but not necessarily the variation index v = w it ^. However, it is exactly an increment in the excess 
parameter which ensures that the bijective map (|7.3ip in the Ramond sector extends to a bijective map 
in the Neveu-Schwarz sector, 



1 „ 1 



C7uC7 ^ U -^fcL- (7-64) 



2 ^2 

"'adm=|f | - ^ 



2 



(7) 

Let us examine Z\ (q) featured in Figure [H From (j7.53p , we have 

2 

,(7)/ s _.JL-L.A„ , [7 



4^^(g) = g-5i+A2,i ^ (7.65) 

2 ^ ^ 
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whereas ()7.54p yields 

2 

= (1 + g + g2 _^ g3 _^ 2/ + 2q^ + + q'^ + q^ + q^ + 

Consistency requires the easily verified relations 

A2,i + 2 = A3,i, A2,i + 5 = A4,i. (7.67) 

With reference to the 21 configurations in Figure O the 14 configurations contributing to ch2^(g) are 
the admissible ones with w = 1. The six configurations with w = 2, which are all admissible, give rise 

(8) 

to ch3^(g), while the single non-admissible configuration has w = 1. According to our combinatorial 
prescription (j7.64p . it contributes to ch^ j^(g) after moving the single occupant in the left column (at level 
3) to the right. Since A2.1 = 1, A3.1 = 3 and A4^i = 6, this is in accordance with the decomposition and 
use of boldface in (j7.46|) . 

We conclude this analysis of the Neveu-Schwarz sector by ignoring the separation into f -sectors. The 
corresponding finitized full partition function is given by 

Z^'^Hq) = Yl ^v'^Hq) = E 2^ch,,i(g), (7.68) 

v=—n r=l 

whose continuum scaling limit yields the full partition function, 

Z(g) = hm Z(^)(g) = Z,{q) = 2Y,rchr,i{q). (7.69) 



In terms of the W-irreducible characters 



Xoiq) = ^;7^(^i'2('z)+r/^('z)) = Yl ^ch,,i(g), Xi{q) = ^^77) (^i'2(g) - ??^(g)) = Y '^chr.iC'Z), 

^^^^ re2N-l rG2N 

(7.70) 

we have 

^ 00 

Xoiq) + xiiq) = ^^^Mq) = Y ^ (^-^i) 

and hence 



Z{q) 

This is in accordance with (|4.28p since ^2(9) 



= 2[xo{q) + Xiiq)]- 
= 2i9i,2(g). 



(7.72) 
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VIII. PERIODIC BOUNDARY CONDITIONS 

In the periodic case, the diagonahsation of was carried out in The relation of the intermediate 
fermions Cj,Cj to the Pauh matrices is the same as in the open case, see (j4.2p . The passage to the 
fermions "Hk^Vk shghtly different, however, and distinguishes the sectors with v even from the sectors 
with V odd. As the notion of sectors is essentially trivial for N odd, our considerations will be split 
according to the parities of N and v. For even, the eigenvectors of have even (resp. odd) fermion 
number in the v even (resp. odd) sector. Because of similarities with the results for open boundary 
conditions, we refer to these sectors as the Ramond {v even) and Neveu-Schwarz (v odd) sectors. 



A. Ramond sector: N even, v even 



For N and v even, in terms of the fermions 1 < k < N, the square of the transfer matrix 

factorizes as [l| 



N_ 

2 



T^ = (^Ak, ^fc = exp|2asingfc ?7fcr/^+i_^| exp|2asing;, qk 

k=l 



{2k - l)7r 
N ■ 



^.11 



The blocks are all four-dimensional and easily diagonalised, and one obtains the relevant spectrum of 



by retaining the eigenvectors with an even fermion number 



identical factors 



rpZ 



k - +1 




implying 








(g) Ak 




(g) Ak 


k=l 




k=l 



•2) 



The first, left, factor in square brackets has a spectrum given by 

|iV±2 

,^ . . , „ . . 2(l-efc-/ifc) 

A; 



f = [y 1 + a^sin^gfc + asinq'fc 



k=l 



Qk 



while the spectrum of the second, right, factor is 



= + o? sin^ gfe + a sin 

k=\ 



2(l-e-fc-/ifc) 



qk 



{2k - l)7r 
N 



{2k - 1)tt 
N 



.3) 



(8.4) 



where the numbers efcj^fcj^fcj/^fc take the values 0,1 as before. We find that the energy levels have the 
following asymptotic form 

c 1 



i>i i>i 

(8.5) 
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with 6j = Ej + fij, 6j = Sj + fij and N' = y- In particular, we see that the conformal spectra of the two 
factors in are identical, and are recognised as two separate Ramond sectors ()7.ip of even system size 
A^', but with vanishing boundary free energy, = (as expected for periodic boundary conditions). 

As is the tensor product of the two factors, it follows that the total energy levels can be written 
in terms of a 'left' and a 'right' contribution, 

= Ef + Ef, (8.6) 

in a way reminiscent of non-chiral conformal representations. The two 'chiral halves' correspond to two 
open Ramond sectors with q replaced by q' = q^-, 

q = e N —^q' = Q -jp- = g n = q . (8.7) 
The variation index v is given as an eigenvalue of the operator 

N 

V = V (-1)^ CjC. = ?7^v+2 V 3N+2 + vIn+2Vn+2 

+ Yl (4^f+fe + ^iv+fe^fc + ^iv+i_fc^JV+l-fe + ^]v+l-fc^f+l-fc)' (^•8) 
and is an even integer in the range 

z;G{-f,-f + 1, ...,f}. (8.9) 

The first two ij^rj terms in (18. Sp are present only for N = 2 mod 4. 

The way the eigenvalues f of V are computed respects the left-right factorisation just given, allowing 
us to write 

V = Vi+Vr, (8.10) 

where 

j j 

Since the variation index v is even, the two chiral parts are not independent, as we have the gluing 
condition 

vi + Vr€ 2Z. (Ramond) (8.12) 

Ignoring this gluing condition, by counting all conformal states, gives rise to a generating function equal 
to the square of that of the open case for even system sizes, namely 0'^{y\q') /rf'{q') where q' = q^. The 
physical partition function in the Ramond sector follows by imposing the Ramond gluing condition, and 
it is therefore obtained by retaining only the terms with even y-powers in the previous function, that is, 
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To determine the contributions from the various f-sectors, Z^{q'), where 



?.14) 



it is convenient to first consider the finitization, 



2LfJ 



i.l5) 



of this expansion following from the separation into chiral halves. To account for all states for a finite 
system size A^, we need to treat the two chiral halves (slightly) differently for N = 2 mod 4, as follows 
from the (conventional) form ()8.2p . We therefore consider that hi > hr, where 



hi=l 



N+2 I 
4 J ' 



i.l6) 



are the maximal heights of the two-column configurations in the left and right chiral halves, respectively. 
We thus have 



where 



Vj^-\-Vr — v 



Vj^-\-Vj-—V 

N-\-6 I I iV+4 



2L 



N+2 



N+2 
- 4 - 



VI 



2LfJ 

N I 



I - 

L 4 



y ^ y ^ Zv;r,f '^^r, 
r=l f=l 



(2[iV±2j+l) (2[f J+1) , 



{q ) ch. 



f,2 



.17) 



Z^,r,f = ^(1 + max{|i;|,r + f} -max{|7;|,|r - f|} 



(8.18) 



This has a stroiig resemblance to the corresponding finitization of the Ramond sector of the CDP model 



on the cylinder [1^, but differs crucially by including also negative variation indices v as opposed to only 
non-negative defect numbers i in the CDP model. A similar difference between the dimer model and 
the CDP model appears in the corresponding Neveu-Schwarz sectors, see (j8.34p . As a consequence, the 
sum of the partition functions in the Ramond and Neveu-Schwarz sectors in the dimer model is modular 
invariant, as discussed below, whereas the similar sum in the CDP model is not modular invariant, as 
discussed in {l3 |. 

In the continuum scaling limit, we see that the partition function in the Ramond sector with variation 
index v is given by 



Z^iq')= E ^^l.,fCh,,2(g')ch,,2(g') 



.19) 
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Ignoring the separation into f-sectors yields the total partition function in the Ramond sector, obtained 
by setting y = 1, 



r=l f=l 

for finite N, and 



Z^{q') = Z\q';l) = Zf;,chr,2{q') ch,,2{q'), (8.21) 



r,reN 

in the continuum scaling limit, where 



Zf;, = E Z^^,^, = 1(1 + i-iy+yr. (8.22) 
In terms of W- irreducible characters (j7.36p , we thus have 

Z''W) = x\{q') + xliq'). (8.23) 

8 8 

This is in accordance with (j8.13p since ^3(5') = '&o,2{q') + '&2,2{q') and 04((7') = '&o,2{q') — '&2,2{q')- 

B. Neveu-Schwarz sector: TV even, v odd 

For even and v odd, in terms of the fermions ?7fc,??|., < A; < — 1, the square of the transfer 
matrix factorizes as 

f -1 

T'^ = {Aq® An)® ^ Ak, ^fc = exp|2asing;t ?7fcr?^_fc| exp|2asin^fc 7?]y_^r?^, qf, = ^. 

k=l 

(8.24) 

The blocks Ak are four-dimensional, while ^0 = I2 and An_ = I2 act in the spaces spanned by {|0), ?7o|0)} 
and {|0), ??!v|0)}, respectively. As for v even, can be written as 



rj-\2 



LfJ L^J 



(X) ^fc 



fc=i ^ fc=i 



(g) Afc . (8.25) 



These two blocks are easily diagonalised, with eigenvalues AJ^^, AJ5^ given by formulae similar to the 
even case, (j8.3p and (j8.4p . For A^ large, we find that the energies are given by 



i>i i>i 
with (5j = Ej + //j, = Ej + /Xj and A^' = -y. Each 'chiral' energy level is doubly degenerate, because of 

^0 = and An_ = I2 appearing in the two blocks. Therefore, the two chiral halves each have exactly the 

2 

same conformal pattern (j7.39p as in the open case of odd system size A^', but with vanishing boundary 
free energy, f^^^ = 0. 
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It follows that the full energies E^^ separate into two chiral contributions, 

^NS^^NS+^NS^ (8.27) 

every such energy level being counted four times. They do not form the relevant spectrum in this sector 
however, as we have to retain the eigenvectors with an odd fermion number, that is with an odd variation 
index. 

The variation index v is given as an eigenvalue of the operator 

TV 

V = V (-1)-^ C]C- = VoVn+ V^n Vo + V^N + vIn 

J = l 

+ Yl ('nkVN+k + ^f+k^k + v\_k^N-k + VN-kVN_k)^ (8-28) 

l<fe<f 

and must be restricted to odd integer values between — ^ and ^. The third and fourth ry^r/ terms are 
present only if = mod 4. 

As in the Ramond sector, the variation index respects the separation into chiral halves and is therefore 
written as 

V = Vl + Vr, VI =Y(.l^j - ^j)^ Vr = Y(fij~^j)^ (^-29) 

j j 

although the attained values are different. In particular, one eigenstate in the block Aq ® An_ contributes 

2 

an additional +1 to v, another one contributes —1, and the last two contribute 0, so that all together, 
they account for a factor {y + y~^ + 2) = (ya + y~'^Y the generating function. In the Neveu-Schwarz 
sector, the gluing condition is 

V = vi+Vr ^27. + !. (Neveu-Schwarz) (8.30) 

Ignoring this gluing condition, by counting all conformal states, gives rise to a generating function equal 
to the square of that of the open case for odd system sizes, namely 02{y\q') / rf' {q') where q' = q^. The 
physical partition function in the Neveu-Schwarz sector follows by imposing the Neveu-Schwarz gluing 
condition, and it is therefore obtained by retaining only the terms with odd y-powers in the previous 
function, that is. 



The contributions from the various v-sectors, Z^^{q'), follow from 

z''H<i';y)= E (8.32) 

whose finitization reads 

2L^J-l 

zNSW(,';y)= E Z:^'^''H<l')y''- (8-33) 

^=_2L^J+l,by2 
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The finitized partition function for fixed v is given by 



{q')zl 



ill 



-\-v.p — v 

Vj^ -\-Vf—V 



-24+^|„;|+l,l 



2LfJ + i 



'NS .J^L^J), ;,{2L^J), , 



(q) 24+'^|«r| + i,l 



2L^J+1 



where 



Z^-rf = max{|u|, r + r} — max{|t;|, \r — r\}. 



In the continuum scahng hmit, we thus have 



r,fgN 



Ignoring the separation into w-sectors yields 



r=l r=l 



and 



where 



r,feN 



VG2Z+1 



In terms of the W- irreducible characters (I7.70p . we thus have 



Z''^{q')=2[Uq') + UQ')?- 
This is in accordance with ()8.3ip since 6i{l\q') = and ^2(9') = 2t9i.2(9')- 



(8.34) 



(8.35) 



(8.36) 



(8.37) 



(8.38) 



(8.39) 



(8.40) 



C. Full partition function and modular invariance for TV even 



The full physical partition function for even is obtained by adding the contributions from the 
Ramond and Neveu-Schwarz sectors, 



Z{q';y) = Z\q'-y)+Z''^{q'-y) 



q' = e —. (8.41) 
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This partition function may be expanded in powers of y to give the partition functions for individual 
sectors. However, this merely yields the decompositions already discussed above. 

We can relax our considerations by not keeping track of the sectors. This corresponds to setting y = 1. 
In this case, we readily recover the torus partition function 

Z{q') = Z{q'; 1) = Z^{q') + Z^^{q') = ^'^'^^ 'J)^ ^'^'^^ (8-42) 



originally obtained by Ferdinand 



331 ]. In terms of W-irreducible characters, this partition function reads 



Z{q') = Xl 1 («') + Xl (q) + 2[Uq') + Xi{q)?- (8.43) 

8 8 



This is modular invariant and also appears as the partition function for symplectic fermions as described 
by the triplet model [10]. 

D. Periodic boundary conditions for TV odd 

The full spectrum for A'^ odd is the same with periodic boundary conditions as it is with open boundary 
conditions, where the open spectrum is given in (j7.39p . The patterns of degeneracies, in particular, are 
thus identical. In the periodic case, each eigenvalue of the form ()7.39p indeed occurs twice; once in the 
even fermionic sector, and once in the odd fermionic sector. It follows that the partial partition functions 
corresponding to an even or odd fermionic number, respectively, are equal, and that the full partition 
function is the same in the odd open and odd periodic cases. We recall that this partition function is 
given by (g^H])- 

As discussed in Section IVj the even and odd subspaces coincide with the two invariant sectors, £^ 
and defined in (j3.16p . Since the two partial partition functions are equal, we need not worry about 
which is which, and accordingly, we write the partition function as 

Z{q) = Z+{q) + Z.iq), (8.44) 

where 

Z+(g) = Z^{q) = ^ = -j^ = Uq) + Uq), Q = e-^. (8.45) 

The corresponding finitizations follow readily from the finitizations in the open case discussed in Sec- 
tion EUll 

IX. CONCLUDING REMARKS 



We have examined the conformal properties of the dimer model based on the working hypothesis 
that the model is described by a conformal field theory with central charge c = —2. Our analysis is 
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significantly finer than previously carried out in the literature, and this improvement was made possible 
by our introduction of a quantum number called the variation index. This quantum number is present 
for finite lattice systems and partitions the set of dimer configurations into sectors. Finitizations of the 
corresponding conformal spectra are described combinatorially in a way reminiscent of similar studies 
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12l | of the critical dense polymer model. Since our results are based on an exact determination of the 
transfer matrix eigenvalues in the various sectors, they are analytical in nature. 

There are still many open problems associated with the dimer model. We hope to return elsewhere 
with (i) a detailed study of finite-size corrections; and (ii) a discussion of the possibility of having non- 
trivial Jordan blocks arising in the continuum scaling limit. Both of these studies are expected to help 
settling the long-standing conundrum with the central charge in the dimer model. We stress in this 
regard that, even though the results of the present work are consistent with a conformal field theory with 
c = —2, they are not definitive. 
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Appendix A: Theta functions 

We use the following convention for Jacobi's theta functions, 

0iiy\q) = --i E (-ir^y''?'^= -iVy 9^ 11(1- (A.l) 

^2 



e2{y\q) = '/'^ = 9^ n (1 - + y'?")(i + y"''?""')' (A-2) 

2 °° 

neZ n=l 

2 °° 

04iy\q) = Yl (-1)" y"?"^ = n (1 - - y9""^)(i - y-'q""-^), (a.4) 

TieZ n=l 

here written in terms of the nomes y and q. At y = 1, they reduce to the standard theta functions, 

ejiq)=ej{l\q), j = 1,2,3,4, (A.5) 



42 



where 9i{q) = and 

^2(9) = ^3(g)=^g#, ^4(g) = ^(-l)-g^. (A.6) 

The related Dedekind eta function is defined by 

00 

V{q)=q^^ll(l-qn- (A.7) 

n=l 

In terms of the generahsed theta functions 
we have 

6^{q) = 2^M, 03(g) = i9o,2(g) + ^92,2(g), 04{q) = ^0,2(9) - ^2,2(9)- (A.9) 

Appendix B: Rewriting of the bulk free energy 

The function arcsinh(a sint) appears in the ground-state energy ()4.12p and in the expression ()4.14p 
for the bulk free energy per site. Representing this function by an integral yields 

TT ZL 1 

T -w-N I ^ f sin t /T-. -, \ 

dt arcsinh(a smtj = j at j dx — j^^^^^^^^= . (^.1) 

V 1 + sin^ t 



.^0 .^0 



Substituting 

n = cost, y = ax, (B.2) 
and performing the integration over n, then gives 



2 



arctan y i 



dt arcsinhfa sin t) = / dy = -(Li2(— ia) — Li2(ia)). (B.3) 

Jo y ^ 

Furthermore, it is recalled that the Legendre chi function X2 is a special case of the Lerch transcendent 

expressible in terms of dilogarithms Li2, 

00 

- ix2(ia) = jH-a', 2, i) = ^(Li2(-ia) - Li^M) = (B.4) 
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